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Abstract 

A light front field theory for finite temperature and density is currently being developed. 
It will be used here to describe the transition region from quark matter to nuclear matter 
relevant in heavy ion collisions and in the early universe. The energy regime addressed is 
extremely challenging, both theoretically and experimentally. This is because of the con- 
finement of quarks, the appearance of bound states and correlations, special relativity, and 
nonlinear phenomena that lead to a change of the vacuum structure of quantum chromody- 
namics. In the region of the phase transition it eventually leads to a change of the relevant 
degrees of freedom. We aim at describing this transition from quarks to hadronic degrees 
of freedom in a unified microscopic approach. 

1 Introduction 

Lattice calculations of quantum chromodynamics (QCD) give firm evidence that nuclear matter 
undergoes a phase transition to a plasma state at a certain temperature T c of about 170 to 180 
MeV. Calculations have been performed at a chemical potential fi = 0. Recent results are for 
staggered fermions JUE1 and renormalization group improved Wilson fermions [3|. The low 
density region reflects, e.g., the scenario during the evolution of the early universe. To achieve 
information from lattice calculations at small // several methods have recently been developed, 
i.e., multiparameter reweighting flUE), Taylor expansion at ji ~ @|7|, imaginary \i ll8ll9lfT0l. 
The region of validity is approximately \i < T ifTTl . Effective approaches to QCD indicate an 
extremely rich phase diagram also for ji > T lfl"2l . Experimentally the QCD phase diagram is 
accessible to heavy ion collisions. In particular relativistic heavy ion collision at SPS/CERN and 
RHIC/BNL explore the region where hadronic degrees of freedom are expected to be dissolved. 
Some results of RHIC are now available that give hints of a non-hadronic state of matter lfT3l . 

On the other hand light-front quantization of QCD can provide a rigorous alternative to lat- 
tice QCD [TT31 . Although the calculational challenge in real QCD of 3+1 dimension seems large 
(as does lattice QCD) it starts also from the fundamental QCD. The light-front quantization of 
QCD has the particular advantage that it is completely formulated in physical degrees of free- 
dom. It has emerged as a promising method for solving problems in the strong coupling regime. 
Light front quantization makes it possible to investigate quantum field theory in a Hamiltonian 
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formulation iflUl . This makes it well suited for its application to systems of finite temperature 
(and density). The relevant field theory has to be quantized on the light front as well, which is 
presently being developed IIB El [IH1 [13 |2Ql ED E3 |23l |M E3 I present the light- 

front field theory at finite temperatures and densities in the next section. For the time being it 
is applied to the Nambu-Jona-Lasinio model (NJL) model ll28l [29 1 that is a powerful tool to 
investigate the non-perturbative region of QCD as it exhibits spontaneous breaking of chiral 
symmetry and the appearance of Goldstone bosons in a transparent way. Finally, going a step 
further I shall give the general in-medium light cone time ordered Green functions that allow us 
to treat quark correlations that lead to hadronization. 



2 Light front thermal field theory 

The four-momentum operator P 11 on the light-front is given by (notation of Ref. El) 

P» = J duo + T + "(x), (1) 

where T^ u (x) denotes the energy momentum tensor defined through the Lagrangian of the 
system and is the quantization surface. The Hamiltonian is given by P~ . To investigate a 
grand canonical ensemble we need the number operator, 

N = J du+j+(x), (2) 

where j u (x) is the conserved current. These are the necessary ingredients to generalize the co- 
variant partition operator at finite temperature ll30ll3Tll32l to the light-front. The grand canoni- 
cal partition operator on the light-front is given by 



Z G = exp jy du+ \-(3 v T +v (x) + aJ+(x)]| , 



(3) 



where a = fi/T, with the Lorentz scalars temperature T and chemical potential ji. The velocity 
of the medium is given by the time-like vector u u u u = 1 OOll . and (3 U = u u /T. We choose 
the medium to be at rest, u u = (u~,u + , u^) = (1, 1, 0, 0). The grand partition operator then 
becomes 

Z G = e- K ' T } K= l -{P- + P + )-fiN (4) 

with P ± and iV defined in ([TJ and ©. The density operator for a grand canonical ensemble ll33l 
l3~4l in equilibrium follows 

PG = (Tr e- K ' T )-' e- K ' T . (5) 

The corresponding Fermi distribution functions of particles f + = f and antiparticles f~ are 
given by 



i n. i 



ex P { 7p ( 2^011 + 2 T ^ ) K 1 



(6) 



and k~ n = (k\ + m 2 ) / k + . This fermionic distribution function (for particles) on the light-front 
has first been given in ifToll . The fermi function for the canonical ensemble can be achieved 
by simply setting fi = 0. This then coincides with the distribution function given recently in 
Ref. 1 21 1 (up to different metric conventions). 



The light-front time-ordered Green function for fermions is 



iGap{x ~y)= 6(x + - y + ) (p a {x)Vp{y)) - %+ - x+) (V p {y)V a (x)) . (7) 

We note here that the light-cone time-ordered Green function differs from the Feynman propa- 
gator Sp in the front form by a contact term 7+ /2k + and therefore coincides with the light-front 
propagator given previously in Ref. [ 35 1 . To evaluate the ensemble average (...)= Tr(pc ■ ■ ■) 
of ©, we utilize the imaginary time formalism ll33ll34ll . We rotate the light-front time of the 
Green function to imaginary value. Hence the A;~-integral is replaced by a sum of light-front 
Matsubara frequencies cu n according to lfT6l . 

1 1 1 _ 1 

— k — > iuJn k + + a = —k„ — > -z, (8) 

where oj n = ttXT, A = 2n + 1 for fermions [A = 2n for bosons]. In the last step we have 
performed an analytic continuation to the complex plane. For noninteracting Dirac fields the 
(analytically continued) imaginary time Green function becomes 

on ~ on 

7fc on + m fl(-fc+) , ^ 7fc on + m g(-fc+) - 

where k = (k + , k±). For equilibrium the imaginary time formalism and the real time formalism 
are linked by the spectral function Il33ll34ll26l . For /i = this propagator coincides with that 
of E6|, but differs from that of ET1B51. 



3 Spontaneous symmetry breaking and restoration 
3.1 NJL model on the light-front 

The Nambu-Jona-Lasinio (NJL) originally suggested in ll2"Sl E91 has been reviewed in Ref. Poll 
as a model of quantum chromo dynamics (QCD), where also a generalization to finite tem- 
perature and finite chemical potential has been discussed. Its generalization to the light-front 
including a proper description of spontaneous symmetry breaking, which is not trivial, has been 
done in Ref. 071 . which we use here. The Lagrangian is given by 

£ = - m )^j + G ($7p) 2 + (^75 1» 2 ) • (10) 

In mean field approximation the gap equation is 

m = m - 2G(W) = m + 2*GA J -0^ TrS F (k), (11) 

where A = NfN c in Hartree and A = NfN c + | in Hartree-Fock approximation, iV c (Nf) is 
the number of colors (flavors). For the isolated case Sp(k) is the Feynman propagator. Taking 
only the lowest order in 1/N C expansion of the 1-body or 2-body operators the light front gap 
equation can be achieved by a k~ integration, where in addition m — > m and G — >• G have 








Figure 1: Effective quark mass as a function 
of temperature and chemical potential. The 
fall-off is related to the vanishing condensate 
(uu), which shows the onset of chiral symme- 
try restoration. Critical temperature at \x = 
is T c ~ 190 MeV. 
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Figure 2: Chiral phase transition as defined 
in The lower part is the chiral broken 

phase, whereas the upper part reflects the re- 
stored phase. 



to be renormalized to accommodate the expansion. For details see ll37l . The propagator to be 
used in (fTTTl is given in ©. The gap equation becomes 

2k+{2n t *rn(T, - /+(*+, k ± ) - /"(*+, k ± )). (12) 

To regularize (fT2l we require k~ a + k + < 2Vt. As a consequence k\ < k + < k\ and 

k\ < 2Qk + - (k + ) 2 -m 2 , (13) 
fc+ 2 = fl=F Vfi 2 -m 2 . (14) 

For the isolated case this regularization is fully equivalent to the Lepage-Brodsky one and the 
three-momentum cut-off. For the in medium case this regularization leads to analytically 
the same expressions as given in [36| for the instantaneous case 11271 . The calculation of the 
pion mass m v , the pion decay constant /„., and the condensate value are also available on the 
light-front 03. 

3.2 Results 

The model parameters are adjusted to the isolated system. We use the Hartree approximation, 
i.e. A = N c Nj = 6. Parameter values are chosen to reproduce the pion mass m n = 140 MeV, 
the decay constant /„. = 93 MeV, and to give a constituent quark mass of m = 336 MeV, i.e. 
G = 5.51 x 10~ 6 MeV, m = 5.67 MeV, and Q = 714 MeV. The parameters are reasonably 
close to the cases used in the review Ref. J36|. The difference between the bare mass m and 
the constituent mass is due to the finite condensate, which is (uu) 1 / 3 = —247 MeV. 

In hot and dense quark matter the surrounding medium leads to a change of the constituent 
quark mass due to the quasiparticle nature of the quark. The constituent mass as solution of (fT2b 
is plotted in Fig. [l] as a function of temperature and chemical potential. The fall-off is related 
to chiral symmetry restoration, which would be complete for mo = 0. It is related to the QCD 



phase transition. For T < 60 MeV the phase transition is first order, which is reflected by the 
steep change of the constituent mass. To keep close contact with the 3M results we have chosen 
for the f2 in-medium regulator mass f2 2 (T, ff) = Ag M + m 2 (T, /i) with A 3 m = 630 MeV fixed 
for all T and [i. 

We define the phase transition to occur at a temperature at which m(T, jj) is half of the 
isolated constituent quark mass [38 1. The phase diagram is shown in Fig. [2 The line indicates 
the phase boundary separating the hadronic phase from the quark gluon plasma phase. Results 
presented in this section are based on an effective interaction in the qq channel. They have to 
be supplemented by the medium dependence of f n and m n that are currently underway. 



4 Few-particle correlations 

We are now interested in the qq channel. Since we are going up to the three-particle system, we 
presently approximate the spin structure. To solve the full three-fermion problem on the light 
front even for the isolated case is quite a challenge. The spin structure is already rather complex, 
see e.g. [39 1. Therefore the elementary spins are averaged Tr7 = (in the medium) and hence, 
for the time being, we are only dealing with bose type particles however subject to Fermi- 
Dirac statistics. Our main focus here is to see how such a three-particle system is dynamically 
influenced by a medium of finite temperature and density; ultimately, how nucleons are formed 
in the hot and dense environment of a plasma of quarks and gluons as the temperature and 
the density becomes smaller and how the relevant degrees of freedom in the Fermi function 
change as the many-particle system undergoes a change to hadronic degrees of freedom. To 
this end we need to formulate suitable few-body equations that describe clusters of quarks in a 
medium. In addition, because of the drastic mass change, see Fig. [H these equations have to 
be relativistic ones. The equations derived here are based on a systematic quantum statistical 
framework formulated on the light front using a cluster expansion for the Green functions. The 
formalism has been given elsewhere [16|. We repeat here the basics to make a connection to 
the previous sections. The light-front time ordered cluster Green function is defined by 

iQafi(x ~y) = 0(x + - V + ) (A a (x)A p (y)) T % + - x + ) (A p (y)A a (x)). (15) 

where all particles A a (x) = A a (x + ,x) = ty ai (x + ,^ty a2 (x + ,x^ty as (x + ,Xj}) ... are all taken 
at the same light front time x + and x = (x + , x±). The upper (lower) sign stands for fermion 
(boson) type clusters. Because of the global light-cone time introduced, the dynamical equation 
for a cluster is equivalent to a Dyson equation with a complicated mass operator that contains 
an instantaneous part and a memory (or retardation) part. For the time being we neglect the 
memory term. This is equivalent to a mean field approximation for clusters leading to Faddeev- 
type three-body equations. 

For a simple zero range interaction the t matrix, Fig. 01 separates and is given by the propa- 
gator t(M 2 ), i.e. 

t(M 2 ) = (iX- 1 - B(M 2 ))~ 1 . (16) 

The expression for B(M 2 ) is represented by the loop diagram of Fig. [Hand, in the rest system 
of the two-body system P M = (M 2 , M 2 , OjJ, given by 

tz>( ]\/f \ 1 f dxd 2 k ± I- f{x, fcj) -/(l -x,kj) 

S(M2) = - (2^ J xjl^) Ml^Mg ' (17) 



Figure 3: Equation for the two-body t-matrix 

with zero range interaction. The crosses refer Figure 4: Loop diagram corresponding to the 
to the Pauli-blocking factor. kernel of the integral equation (El). 



where M| = (k\ + m 2 )/x(l — x) and f = f~ given in © with x = k + / P 2 + where P 2 = 
k 1 + k 2 . For a fermi system there are two important effects occurring due to the blocking factors 
of (flTt . One is the dissociation limit (Mott effect) where M 2 (T d , fj, d ) = 2m(T d , n d ). Above a 
certain temperature and density no bound states can be formed. The second effect is related to 
the appearance of a bose pole in the t matrix. This happens for M 2 (T C , /i c ) = 2/i c and defines 
the critical temperature below which the system becomes unstable and forms a new vacuum 
consisting of Cooper pairs or a condensate. This is related to superconductivity or superfluidity. 
In this case for M| — > Mq 2 , we get 



/(•-a-i) _ . = /(i-*,* 



2 



1 



(18) 



i.e. both nominator and denominator of (fTTt are zero. 

The three-particle case is driven by the Fadeev-type in medium equation 

r / — * \ i / dxd?k ± l-f(x,k 2 ± ) -f(l-x-y,(k + g)l) - 

nv,1±) = (2^)3 W 7 ^i-y-x) M|3m| r ^'^' 

(19) 

where we have introduced vertex functions T and t(M 2 ) given before, and an invariant cut- 
off Mf < A 2 . Here the mass of the virtual three-particle state (in the rest system P M = 

(M 3 ,M 3 ,0 ± )is 

Mqo = S + ^ + a + m2 + g + g 1 + m2 , (20) 
x y 1 — x — y 

which is the sum of the on-shell minus-components of the three particles. The nucleon scale 
is introduced by setting M 3 = 938 MeV. The isolated quark mass used in these calculations is 
m = 386 MeV. 

Fig.|5]shows a shaded area that reflect the region where the transition from baryons to quarks 
(or quark diquarks) occur. The area is defined by use of different regularization masses. The 
chiral phase transition given before is indicated by the solid line. Fig. |6] shows the possible 
transition of quark matter to a superconducting phase. 



5 Conclusion and Outlook 

We have given a relativistic formulation of field theory at finite temperatures and densities uti- 
lizing the light front form. The proper partition operator (and the statistical operator) have been 
given for the grand canonical ensemble. The special case of a canonical ensemble is given for 
/i = 0. The resulting Fermi function depends on transverse and also on the k + momentum com- 
ponents. The k + components emerge in a natural way in a covariant approach. As an application 
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Figure 5 : Nucleon dissociation region (shaded 
area due to different cut-offs). Solid line chiral 
phase transition of Fig. 12 



Figure 6: Phase diagram supplemented with 
critical temperature for color superconductiv- 
ity. Different cut-offs: A = 4m (dot), A = 6m 
(dash), A = 8m (dash-dot). 



we have revisited the NJL low energy model of QCD. We reproduce the phenomenology of the 
NJL model, in particular the gap-equation and the chiral phase transition. We have further given 
consistent relativistic three-quark equations valid in a dense medium of finite temperature. We 
find that the dissociation transition and the critical temperature for the color superconductivity 
agree qualitatively with results expected from other sources. However, the latter results are by 
no means final. We have shown that it is possible to write down meaningful consistent equations 
to solve the relativistic in-medium problem on the light front. The next steps would be to use 
the NJL model all through to give a consistent picture for the qq and the qq channel. Further 
insight into this just emerging possibilities of treating relativistic many-particle systems on the 
light front might be provided by other theories, like 1+1 QCD, the Yukawa model, and finally 
real QCD. 
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